Theory of the Thermal Hall Effect in Quantum Magnets 
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We present a theory of the thermal Hall effect in insulating quantum magnets, where the heat 
current is totally carried by charge-neutral objects such as magnons and spinons. Two distinct types 
of thermal Hall responses are identified. For ordered magnets, the intrinsic thermal Hall effect for 
magnons arises when certain conditions are satisfied for the lattice geometry and the underlying 
magnetic order. The other type is allowed in a spin liquid which is a novel quantum state since 
there is no order even at zero temperature. For this case, the deconfined spinons contribute to the 
thermal Hall response due to Lorentz force. These results offer a clear experimental method to prove 
the existence of the deconfined spinons via a thermal transport phenomenon. 
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The ground state and low energy excitations of corre- 
lated electronic systems are the subject of recent inten- 
sive interests, and especially the possible quantum liquid 
states are the focus both theoretically and experimen- 
tally [lj-Q- F° r the quantum magnets, magnetic suscep- 
tibility, neutron scattering and specific heat are the ex- 
perimental tools to study this issue. In the conducting 
systems, on the other hand, charge transport properties 
also offer important clues to the novel electronic states 
such as the non-Fermi liquid or the quantum Hall liq- 
uid. Therefore, a natural question is whether there are 
any transport properties in insulating quantum magnets 
which provide insight into the ground state. To answer 
this question, we study in this Letter the thermal Hall 
effect theoretically and find several different mechanisms 
leading to the classification of the quantum magnets. 

For a finite Hall response, time-reversal symmetry 
must be broken due to magnetic field and/or magnetic 
ordering. The Hall effect in itinerant magnets, where 
the spin structure and conduction electron motion are 
coupled, has been studied extensively. In this case, in 
addition to the usual Lorentz force, the scalar (spin) chi- 
rality defined for three spins as S% ■ (Sj x Sk) plays an 
important role (fHH . The scalar chirality acts as a ficti- 
tious magnetic flux for the conduction electrons and gives 
rise to a non-trivial topology of the Bloch wave functions, 
leading to the Hall effect. It is natural to expect that a 
similar effect occurs even in the localized spin systems 
for e.g. the spin current Q. Another important tool 
to detect charge-neutral modes is the thermal transport 
measurement. In low-dimensional magnets, the ballistic 
thermal transport property was predicted from the in- 
tegrability of the one-dimensional Heisenberg model 11(1 
and has been experimentally observed in Sr2Cu03 [111 ]. 
In k-(ET) 2 Cu2(CN) 3 , one of possible candidates for two- 
dimensional quantum spin liquids [l[ , the thermal trans- 



port measurement was used as a probe to unveil the na- 
ture of low-energy spin excitations [12| ]. The measure- 
ments have been limited to the longitudinal thermal con- 
ductivity so far. In this paper we predict a non-zero ther- 
mal Hall conductivity, i.e., the Righi-Leduc effect which 
will provide important information as described below. 

First, we need to consider the influence of the external 
magnetic field on localized spin systems. In addition to 
the Zeeman coupling, we have the ring exchange process 
leading to the coupling between the scalar chirality and 
external magnetic fields. This coupling is derived from 
the t/U expansion for the Hubbard model at half filling 
with on-site Coulomb interaction U and complex hopping 
Til 14 1 and its explicit form is given by 
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sin* Si • (Sj x Sk), 



(1) 



where $ is the magnetic flux through the triangle formed 
by the sites i, j, and k in a counterclockwise way. Since 
the coefficient is proportional to t 3 /U 2 , it is expected to 
be small. In the vicinity of the Mott transition, however, 
this coupling is not negligible. We first examine the effect 
of .ffring within the spin-wave approximation. Then we 
find that if the lattice geometry and the magnetic order 
satisfy certain conditions, the magnons can experience 
the fictitious magnetic field and there occurs the intrin- 
sic thermal Hall effect, i.e., the thermal Hall conductivity 
K xy due to the anomalous velocity of the magnons. In this 
case, k xv is independent of the lifetime of magnons 
whereas the longitudinal one k xx depends on r 
It can be regarded as a bosonic analogue of the quan- 
tum Hall effect with zero net flux [171 . We also derive 



a TKNN-type formula [18( of the thermal Hall conduc- 
tivity for a general free-bosonic Hamiltonian. It should 
be possible to apply this formula to the recently found 
phonon thermal Hall effect fl9j | . Finally, we consider the 
effect of -ffring in quantum spin liquids. Since there is no 
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FIG. 1: (a) A portion of the triangular lattice with ferromag- 
netic order. The blue arrows indicate the directions of spins. 
The counterclockwise rotations (red arrows) indicate the or- 
der of sites in H r i n g- We use this convention throughout this 
Letter, (b) Antiferromagnetic model on a square lattice. The 
interaction along the diagonal edge is smaller than that along 
the vertical or horizontal edge, (c) A portion of the triangular 
lattice with the 120° order. 



sharing a link occurs for the collinear antiferromagnetic 
configuration. Finally, for noncollinear spin structures, 
we considered the 120° magnetic order on a triangular 
lattice (Fig. (He)) and conclude that the cancellation 
again occurs since ordered components of S{ and that of 
Si are the same as shown in Fig. QJc). 

Intrinsic thermal Hall effect in the spin-wave approx- 
imation. - Once we understand the principles of the 
cancellation, it is rather easy to find an example where 
it does not occur, namely, when the link is shared by in- 
equivalent cells. An example is the ferromagnetic model 
on the kagome lattice. In this case, the spin wave Hamil- 
tonian is influenced by the magnetic flux $. We will 
develop below a theoretical formalism to calculate the 
thermal Hall conductivity in terms of the Kubo formula. 
For this purpose, we consider a general Hamiltonian for 
noninteracting bosons which can be regarded as a spin- 
wave Hamiltonian within quadratic order in 6S: 



magnetic order in such a system, it has been proposed 
that deconfined fermionic spinons exist. In contrast to 
the magnons, the spinons, which are the gauge depen- 
dent object, can feel the vector potential A just as in 
the case of electrons, leading to the Landau level forma- 
tion [T3|. We propose a novel way to detect the spinon 
deconfinement via the thermal Hall effect measurement in 
a candidate of quantum spin liquid, k-(ET)2Cu2(CN)3. 

NO-GO theorem for the coupling to magnetic flux. - 
Let us first consider the spin-wave expansion of Eq. ([T]) 
to find a system in which the intrinsic thermal Hall ef- 
fect occurs. We consider the collinear ground state spin 
configurations. The fluctuation of the scalar chirality up 
to the second order in SS = S — (S) is written as 

(Si)-(5Sj x SS^+iSj) ■ (SS k x SSi)+(S k ) -(SS t x SS,), (2) 

where (S a ) denotes the a-component (a = x, y, or z) of 
the ordered moment. Note here that the linear order 
terms in 5S vanish since (Si) x (Sj) = 0. As an exam- 
ple, we consider the ferromagnetic Heisenberg model on 
a triangular lattice shown in Fig. QJa) with an ordered 
moment So along y. In this case, the quadratic terms 
in 5S from Eq. (TTJ) always cancel. To explain this, let 
us focus on the edge (jk). From the upper triangle, this 
edge gives So(SSj x SSk) y - On the other hand, from the 
lower triangle, it gives So{SSk x SSj) v which cancels out 
the former one. Since such a cancellation occurs on any 
edge, -f/ring in Eq. (JTJ does not contribute to the spin- 
wave Hamiltonian to quadratic order. This observation 
leads us to conclude that such a cancellation occurs for 
any ferromagnetic model where each edge is shared by 
the equivalent cells such as plaquettes and triangles. A 
similar cancellation occurs for certain antiferromagnetic 
systems. An example is shown in Fig. Gib). In this ex- 
ample, there are several different types of ring exchange 
processes, but again the cancellation between the cells 



H — h(Rj a ), h(Rj a ) — 



-Vt? &t ? b n 



h.c, 



where annihilates a boson at the lattice point of the 

a-th site in the j-th unit cell and 5 a are vectors connect- 
ing Rj a and its neighboring sites. The hopping tg is in 
general complex. In momentum space, the Hamiltonian 
is written as H — J^kbaifyWa^ifybpik) , where b a (k) 
is a Fourier transform of b^ and repeated indices are 

summed over. Using the local energy operator h(Rj a ), 
the average energy current density is defined by 



j E = i[H,J2R ja h(R ja )]/V, 



(3) 



J ot 



where V is the total volume [20|. Using the 
Fourier transform of J\ h(Rj a ) defined by h(q) = 
a e lq ' Rja h(Rj a ), the energy current density is rewrit- 
ten as je — dq[h(0), h(q)]\f =0 /V ■ Using this fact, the 
following convenient expression for je is obtained: 



(4) 



where the differential operator acts only on W(k) 2 . We 
introduce the spin wave basis \u a (k)) which diagonalizes 
T-L(k) with eigenvalues u> a (k). It is important to note that 
even in this basis je is not diagonal. In addition to the 



expected diagonal term 



duj a 



t (k), there are off-diagonal 



terms which can be thougiit of as arising from anomalous 
velocities. As we see below, these terms are responsible 
for k xv , just as in the case of the intrinsic anomalous Hall 
effect in metals 2l|. 

Starting from the Kubo formula, the following expres- 
sion analogous to the TKNN formula [l8[ can be obtained 
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FIG. 2: (a) The unit cell of the kagome lattice (the area 
enclosed by the dotted line). The vectors Si (i — 1, 2, and 3) 
correspond to green, red, and blue arrows, respectively. The 
arrows on the edges indicate the sign of the phase factor e 1 '*'' 3 
(see main text). The fictitious fluxes through triangle and 
hexagon are (j> and —20, respectively, (b) Dispersions of bands 
in the case of <j> = 7r /3 with JS = 1. The right figure shows 
the enlargement of dispersion around uj = where the lowest 
band has a quadratic dispersion, i.e., wi(fe) ~ JS(k^. + ky). 



for the thermal Hall conductivity n xy : 

d 2 k 



K = - 



— ImV 
2T ^ 



a(k) 



bz(2^) 2 

x (a fcx u a (k) I (H (k) + uj a (k) f | d ky u a (k) ) , (5) 

for the non-interacting spin waves (free bosons) where 
the integral is over the Brillouin zone (BZ), and n a {k) = 
^ e l3ui a {k) _ j g ^ e g ose distribution function. We can 
show that the integrand in Eq. ([5]) is independent of a k- 
dependent phase change of the spin wave basis functions 
\u a (k)) in a similar way as for the orbital magnetiza- 
tion |22j |. 

Thermal Hall effect in kagome ferromagnet. - We 
now apply the above formula to the ferromagnetic model 
on the kagome lattice. The Hamiltonian is given by H = 



A,v 



h>A + h-u with 



h 



a/v — ~J{Si ' Sj + Sj ■ Sk + Sk ■ Si) 



jSi-iSjxSk), 



where —J is the ferromagnetic exchange coupling, K 
is proportional to sin$ according to Eq. (J]), and the 
sum is taken over all the triangles in the kagome lat- 
tice (see Fig. EJa)). Again, note that the sites i, 
j, and k form a triangle (A or y) in a counterclock- 
wise way. Using the Holstein-Primakoff transformation 
[5+ = (25 - njfPb^Sj = 6} (25 - n,) 1/2 , 5| = 5 - n, 
with rij — bjbj], we obtain the spin- wave Hamiltonian as 



H sw = 4JSj2 n j - S\/J 2 + K 2 J2 (e~^ /8 6j6* + h.c), 



where tan(</>/3) = K/J and the sum is taken over all 
the nearest neighbor bonds. This Hamiltonian con- 
tains a non-trivial phase factor and is quite similar to 
the electronic model considered by Ohgushi et al. 
The Fourier transform of the Hamiltonian is given by 
%{k) = US - 2 J5(cos(0/3))~ 1 A(fc, <j>) with 

cosfcie-'^ 3 cosfc 3 e^/ 3 

A(k,<f>) = ( cosfcie^/ 3 cosfc 2 e- ,; ^ 3 

cosfc 3 e-^/ 3 cos/c 2 e l0/3 

where ki = k ■ a, with a% = (— 1/2, — V3/2), a,2 = (1,0), 
and 03 = (— l/2,v / 3/2) as shown in Fig. [2a). If 
the magnetic field is absent, i.e., a = (f> = 0, the up- 
per band is dispersionless and the other two bands are 
particle-hole symmetric around u — 3J5. The disper- 
sions of three bands (0 < 0Ji(k) < w 2 (fc) < u}^{k)) 
for <fi — ir/3 are shown in Fig. HJb). In the limit 
of low temperature and weak magnetic field, the dom- 
inant contribution to the integral in Eq. (|S|) comes from 
a = 1 (lowest band) and small \k\ due to the Bose fac- 
tor n a (k). By an explicit calculation [23j |. we find that 

(0fc s Ui(fe)|0fc B i(/)\k\ 2 /(27V3) around fc = and 

obtain 



(6J5) 2 



dk 



k /-# 2 \ 



TT(j} 



2T J 2^e^ fc2 -l U7V3/ 36^3 



T, 
(6) 

where we have replaced the integration over the BZ with 
that over all k. In this way, a non-zero thermal Hall con- 
ductivity is indeed realized by the coupling between the 
scalar chirality and the magnetic field in the ferromag- 
netic kagome lattice. The obtained linear dependence of 
K xy on T is due to the dimensionality (D = 2) and the 
quadratic dispersion around cj(k) = 0. 

Thermal Hall effect in quantum spin liquids. - As 
discussed above, the spin Hamiltonian contains the mag- 
netic flux and the spin waves or magnons are influ- 
enced by the magnetic field only through the off-diagonal 
matrix elements of the thermal currents, correspond- 
ing to the intrinsic anomalous Hall effect. This is in 
sharp contrast to the case of electrons which is cou- 
pled to the vector potential A, and the usual Hall ef- 
fect due to the Lorentz force occurs there. In this re- 
spect, it is interesting to note that the spin operator Si 
can be represented by the fermion operators f\ a and f %o 

(called spinons with spin 1/2) as 5, = J2a a' fur&fiv'/^ 
(<? = (a x , o~y, o~ z ): Pauli matrices, a, a' =t,4) with the 
constraint ^2 a f\ a fi a — 1- In this representation, the ex- 
change interaction J Si ■ Sj can be written as —JxljXij/2 

with Xij = Yjafitrfjv Tnis Xrj is called the order pa- 
rameter of the resonating valence bond (RVB) which de- 
scribes the singlet formation between the two spins Si 
and Sj . In the mean field approximation for Xij i the free 
fermion model for fj a , fi a emerges [24j ]. The phase a^- 
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of the order parameter (xi 



\(Xij)\< 



and the La- 



grange multiplier a to impose the constraint above con- 
stitutes the gauge field, which is coupled to the spmons. 
In the confining phase of this gauge field, two spinons 
are bound to form a magnon. On the other hand, in 
the deconfining phase, the spinons behave as nearly free 
quasi-particles. The latter case is realized in some of 
the quantum spin liquid states [24]. Similar state has 
been obtained also for the Hubbard model 25 1, which 



contains the gapped charge excitations. This charge de- 
grees of freedom is represented by the U(l) phase factor 
e l9i , i.e., the electron operator a a is decomposed into 
the product fi a G~ l6i , which is coupled to a,^ — A^j where 
Aij is the vector potentialfPeierls phase) corresponding 



to the magnetic flux $ 
A = (l/5)/*E P (^' 



F ) 2 



Then, the Maxwell term 
(g: coupling constant, 
d„a^, F^ v = d^Ay - d v A^ ) is generated 
by integrating over the charge degrees of freedom. To 
summarize, the spinons are described by the Lagrangian 



(7) 



Following the previous works [LJ, |25j , we take the spinon 
metal with a Fermi surface as a candidate for the 2D 
quantum spin liquid realized in K-(ET)9Cu9fCN)?i[l2l|. 
In a magnetic field F xy = B z , the average of the gauge 
flux (J^xy) = cF xy is induced with c a constant of the 
order of unity because of the coupling between T xy and 
F xy in C g [lj, [23]. Therefore, the spinons are subject 
to the effective magnetic field (F xy ) and to the Lorentz 
force. 

Let us first estimate the spinon lifetime r from the 
recent thermal transport measurements in that material. 
The longitudinal thermal conductivity is obtained from 
the Wiedemann- Franz law by assuming a Fermi liquid of 
spinons: 



1 



xx _ o 71 " 2 f £ F \ k l T 
~ l 3 \h T ) h 'd' 



(8) 



where e-p is the Fermi energy and d ~ 16A is the inter- 
layer distance. After a subtraction of the phonon COn- 



SI 



tribution, 
T = 0.3 

£F = J ~ 

we examine k 



is estimated to be ^0.02 WK m 



at 



K [12j. We obtain eft/H — 56.5 and with 
250 K, we estimate r ~ 1.72 xLO" 12 s. Next 

the gauge 



As has been shown in [lJl 



flux for spinons is comparable to the applied magnetic 
flux and hence K x £ ~ (uj c t)k xx , where w c = eB/m c is 
the cyclotron frequency with the effective mass of spinon 
m c . Estimating m c ~ h 2 /(Ja 2 ) with assuming the lat- 
tice spacing a ~ 10A, we obtain uj c t ~ 0.086 x B with 
B being measured in Tesla. Therefore, the thermal Hall 
angle K^I K sp ~ w c t becomes of the order of 0.1, which 
is easily measurable, with a weak magnetic field B ~ IT 



such that the spin-liquid ground state is not disturbed. 
Also note that compared with the intrinsic thermal Hall 
effect discussed above, the magnitude of this Lorentz- 
force driven thermal Hall conductivity is much larger by 
the factor of ~ {syt/K) 2 . Therefore, the observation of 
the thermal Hall effect is a clear signature of such de- 
confined spinons in the spin liquid, and experiments on 
k-(ET) 2 Cu 2 (CN) 3 is highly desirable. 

Another important difference between the spinon con- 
tribution and the intrinsic term is that the spinons are 
diffusive and see the field A. Thus in a small sample 
one can expect mesoscopic effects such as universal con- 
ductance fluctuations of the thermal conductivity as a 
function of B. Using the Wiedemann- Franz law, we ex- 
pect the relative fluctuation in k xx and K xy to be of order 
h/ (eft) for each coherent volume with dimension \J Ui n 
where £ = v^t and £i n = vpTi n and 7j„ is some inelastic 
scattering time much longer than r at low temperatures. 
The fluctuation is reduced by y/~N if the sample contains 
N coherent volumes. We estimate the elastic mean free 
path £ to be 400 A, so that at low temperatures this effect 
may be observable in micron-scale samples. 

In conclusion, we have studied theoretically the ther- 
mal Hall effect in the quantum spin systems induced by 
the external magnetic field. There are three cases, i.e., 
(i) no thermal Hall effect, (ii) intrinsic thermal Hall ef- 
fect by the magnons, and (iii) large thermal Hall effect 
due to the Lorentz force, (i) corresponds to the most 
of the conventional (anti)fcrromagnets on triangular and 
square and cubic lattices, while (ii) to the magnets on a 
particular lattice structure such as kagome, and (iii) to 
the spin liquid with deconfined spinons. Therefore, the 
thermal Hall effect offers a unique experimental method 
to gain an important insight on the ground state/low en- 
ergy excitations of the quantum magnets. 
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Supplementary Items for Theory of Thermal Hall Effect in Quantum Magnets 

Off-diagonal terms in energy current density 

In this section, we give an explicit expression for the energy current density je in the basis where the Hamiltonian 
is diagonal. As we have explained in the main text, je is not diagonal in this basis. We start with the non-interacting 
Hamiltonian in the Fourier space: 

h=J2 bUk)n a Ak)b p (k), (9) 

fc,a,/3 

where b a (k) denotes the Fourier transform of and repeated indices are summed over. The single-particle 

Hamiltonian H(k) can be diagonalized by the unitary matrix g(k) as g> (k)H(k)g(k) = Sl(fc), where f2(fc) = 
diag(tJi(fc), Lj n {k)). Using g(k) we define a new basis of bosons as 

3 

In this basis, H is written as H = w a (fc)7^(fe)7 Q (fc). 

Let us now see the off-diagonal structure of je in the basis of 7 a (fc). This is clarified by introducing the Berry 
connection (or Maurer-Cartan 1-form) defined by Q(k) — g^ (k)d^g(k). Using 6, the energy current density (Eq. (4) 
in the main text) is written as 



k 



up 



Here we have used the relation (d?g* '(k))g(k) = —g'(k)(drg(k)). The first term in the bracket is a diagonal one 

corresponding to 9u ^ k ^ uj a (k). On the other hand, the second term represented by the commutation relation gives 
the off-diagonal elements of the energy current density. Similarly to the case of the intrinsic anomalous Hall effect, 
this off-diagonal term can be regarded as arising from anomalous velocities. The thermal conductivity tensor can be 
obtained by substituting the above expression for je into the Kubo formula: 



V 



-I dt I d\(f E (-i\)j E (t)) th , (12) 



o 



where j3 — l/T is the inverse temperature and (...)th denotes thermal average. Along the same lines as the derivation 
of the TKNN formula PI, one can obtain Eq. ©. 
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A detailed derivation of the thermal Hall conductivity in the Kagome ferromagnet 

Eigenvalue problem 

In the main text, we have studied the thermal Hall effect in the kagome ferromagnet. To obtain the spin- wave 
spectrum and eigen-modes, we need to solve the eigenvalue problem associated with the following matrix: 

1 cos(fc-ai)e- 1 ^ 3 cos(fc • a 3 )e 4 */ 3 

A(k,(f>)= cos(fc-ai)e^/ 3 cos(fc ■ a 2 ) e -^/ 3 |. (13) 

y cos(fc ■ a 3 )e-^/ 3 cos(fc ■ a 2 )e^/ 3 

Before study this matrix, we consider a more general matrix 

/ c 

M=\ a* b I (14) 




and its secular equation given by 

A 3 = (H 2 + \b\ 2 + \c\ 2 )X + (abc + a*b*c*). (15) 
Now we apply Viete's solution of cubic equation. To do so, we compare Eq. (fT5|) with the following identity: 



cos3 (!) = i cos (!) + \ cos e - (i6) 

Then, we obtain one of the solutions Ai = 2Acos(6 /3) with cos (9 = B/{2A) (0 < 8 < tt), where A and B are defined 
through the following relations: 

3A 2 = \a\ 2 + \b\ 2 + \c\ 2 , A 2 B = abc + a*b*c*. (17) 

The other two solutions can also be written in terms 9 as A 2 or A3 = 2v4cos[(6* ± 2ir)/3). The normalized eigenvector 
corresponding to A„ (n = 1, 2, 3) is explicitly given by W~ 1 (X 2 i - \b\ 2 , X n a* + be, X n c + a*b*) T where W 2 = 3[2A 2 A 2 + 
{A 2 -X 2 n )\b\ 2 + X n A 2 B]. 

Let us now apply the above general result to our specific problem of A(fc, </>). We identify 3^4 2 with 1 + f(k), and 
A 2 B with f(k) cos(0), where f(k) = 2 Il?=i c °s(£ • Sj). Then, the eigenvalues of A(k, </>) are given by 



X n (k, 0) = 2 J i±M cos [ ^ ] , (n = 1, 2, 3) (18) 



with 



cosfli(fc) = W (f( k ) co l<t>) 2 e 2 (k)=6 1 (k)+2n, 6 3 (k) = O^k) - 2tt. (19) 

V 4 (i + /(fc)) 3 

It is useful to note that 1 + f(k) — 53 7 *=i cos 2 (fc • The normalized eigenvector corresponding to X n (k, </>) is written 

X n (k, <t>) 2 - cos 2 k 2 

\u n {k)) = U I e^/ 3 A„(fc,<?!))cosfci + e~ 2#/3 cos fc 2 cos fc 3 |, (20) 

W n (k,(f)) y e -io4/3 An (^0) cosA . 3 + e 2#/3 cosfclCosfc2 

where W n (k,(f)) 2 = 2(1 + f(k))X n {k) 2 + (1 + f(k) - 3A„(fc) 2 ) cos 2 fc 2 + 3A„(fc)/(fc) cos 4>. As for the ferromagnetic 
Kagome model, the Hamiltonian is defined by W(k) — 4JS — 2 JS(cos(4>/3))~ 1 A(k, 0) and hence the eigen-energies 
are given by 



< , . l + /(fc)cos(g»(fc)/3) , 

Wn (fc)=4j5|i-y— — ^1- ( 21 ) 
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Derivation of Eq. (6) 



We now derive Eq. (6) in the main text from the general formula for n xy (see Eq. (5)). In the limit of low 
temperature and weak magnetic field, the contribution from the lowest band (a = 1) dominates in RHS of Eq. (5) 
due to the Bose factor n a (k). Therefore, we approximate Eq. (5) as 

d 2 k 



2T 



Tm 



BZ 



(2tt)5 



ni(fc) (dk x ui(k)\(H(k) + u>i(k)Y\dk y ui{k)) . 



(22) 



Inserting the resolution of the identity at each fc, i.e., 1 = X^=i \ u p{k)){ u p(k)\, we can rewrite the above expression 



as 



— V 

2T ^ 



(9=2,3 



BZ 



d 2 k 



ni(k)(u)i{k) + w /3 (fc)) 2 Im (<9 fcx wi(fc) up(k)) (up(k) d ky ui{k) 



(23) 



Note that we have omitted /3 = 1 from the summation above since it can be shown that 
{dk w u\{k)\ui(k)){u\{k)\dk y ui{k)) is real using the fact \ui(k)} is normalized, i.e., {d^u\{k)\ux{k)) + {u\{k)\d^ui(k)) — 
0. In the low temperature limit, the dominant contribution to the integral in Eq. ([23]) comes from small |fc| due to 
the Bose factor. In the vicinity of k = 0, one can expand cos#i(fc) in Eq. (fT9|) and Ai(fc, </>) in Eq. (fT8| with respect 
to k as 



cos 6»i (/c) ~ cos 0(1 + 0(fc 4 )), Ai(fc,( 



2cos(0/3) l^~\^y i J +0{k i )^ 



(24) 



Therefore, in the vicinity of k = 0, 0Ji(k) ~ JS'lfcl 2 . Furthermore, using Eq. (fT9|) and (p4|) . we find wi(fc) + W2(fc) ~ 
2 J5(3 - \/3 tan(0/3) + 0(fc 2 )) and wi(fc) + oj 3 (k) ~ 2 J5(3 + v / 3tan(0/3) + 0{k 2 )). Then, we approximate uj 1 + uj 2 
and wi + LU3 by 6 JS and rewrite Eq. (|23[) as 

d 2 fc 



1 

2T 



BZ 



(2tt) : 



(k)(6JS) 2 J2 Im ( d k*Mk) up{k))(up{k) d kyUl {k) 



(3=2,3 



(6JS) 2 
2T 

(6JS) 2 
2T 



d 2 fc 



■ni( 



(fc)^Im (dk^u^k) up{k)\ (up(k) d ky ui(k) 



(8=1 



BZ 



■ni(fc) Im K8 fc(D ui(fc) d ky ui{k) 



(25) 



Here, we have again used the fact that (9fc x «i(A;)|Mi(fc))(ui(fc)|9fe !J Mi(fc)) is real. Note that if we take into account the 
(ft- and fc-dependencies of uj\ + Ui (i = 2, 3), they just give higher order terms in <fi and T. (As we will see in the next 
paragraph, hxi[(dk m ui(k)\d ky ui(k))} is already first order in 0. ) 

For our purpose to obtain the thermal Hall conductivity, we need to calculate the imaginary part of 
(dk m ui(k)\dk y ui(k)} . Now we consider the perturbative expansion of \ui(k)) with respect to <fi an d keep up to the 
first order terms in 4>. The zeroth and the first order wavefunctions are given by 

/ A 1 (fc,0) 2 -cos 2 /c 2 

\\(k, 0) cos/ci + cosfc2 cos &3 I (26) 
\ Ai (k, 0) cos k 3 + cos k\ cos k 2 



1 

Wi(A?,0) 



and 



| U (D 0^ = K^/ 3 ) I Al o) cos fci - 2 cos k 2 cos k 3 
Wi(fc,0) ^ _Ai(fc,0)cosfc 3 + 2cosfc lC osfc 2 

respectively. Then, from the Taylor expansion around k — 0, we obtain 

lm[<0 fc .u 1 (fc)|0 fc ,u 1 (fc))] ~ M^k."! 05 ®^,"^*)) + {d k ^\k)\d kv uf\k))\ 



(27) 



(28) 
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Therefore, for the first order term in cj>, the Berry curvature around the point k = is quadratic in Substituting 
Eq. (|28p into Eq. (1251) . we obtain Eq. (6) shown in the mai n text. 



